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ABSTRACT 



j It is common to estimate a distribution by means of a step function. 
Such estimates can be made continuous by connecting the points of the 
steps with straight line segments. In this paper the best estimator 
of this class is found for data which is uniformly distributed using 
minimum risk, then this risk is compared with those of the sample 
distribution function and the Pyke estimator. 
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I ^ m'CRODUCTION 



In reference > Prof. Read studies the properties of some estimators 
for a continuous distribution. In particvilar , he computes the risk function 
(xising the assTomption of the vmiform distribution) of the estimator that 
connects the points C the r smallest obsesrvation , 

vTith straight line segments . This estimator is shorn to be better than 
the aample distribution function in an asymptotic senae . It is called 
the Pyke estimator , 

It is natural to ask the question , "■ How well can a distribution be 
estimated by a f\anction that connect the points with straight line 

segments?" The purpose of this paper is to find the sequence that 

minimizes the integrated risk function (using the uniform distribution) 
and to compare that risk function with those of the sample distribution 
function and the Pyke estimator. 
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IIJ)ETERMIMATIOU OF THE OPTIHDM C0EFPICIE1IT3 



Let ^ increasing sequence with. 0^ = 0 and 

= 1 .A continuous function estimator for P(x) can be defined by 



(2.1) H(x) = C 



x-X 






I ,-x ■ 

r+1 r 



for ; X X . r = 0,..,n 

’ r r-f1 * ' 



where order statistics of a random sample 

from an absolutely continuous distribution function P(x), It is convenient 
to assume that the popiilation sampled is bounded and contained in a 

X . = b and the 

n+1 

risk function (R) by 



finite interval [a,.] •Thus, we can define 



^0 = 



-J’ r 

(2.2) R = E J P(x) - H(x) I dP(x) 

a *- “* 

For convenience the interval ra,b] can be taken to •'fh® 

estimator E(x) is defined piecewise according to v/hich of the random 

intervals ) contains x . Assiaming the data come from a uniform 

distribution The joint p.d.f. of any two order statistics , say X and 

r 

» is as easily espressed . Using X^ = u and X^^.j = v the joint 
p.d.f. is 
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n! . n T> 1 

(2.3) fy Y (u,v) - “U ""aI-v) for 0 ^ r ^ n 

r’ T4-1 (r-1 ) ! (n-r-1 ) ! 



and 0 ^ u ^ V 1 






(2.4) y (0,v) - n(1 


xn-1 

- v) 


for r = 0 and 


(2.5) X " nu^"^ 

n* n-hl 


for r = n and 


Let define 






(2.6) 


It 

o 

+ 


- C 

r 


So that 






(2.7) 


II 


A. 

J 


and for simplicity f^ ^ 


(u,v) = fj.(u,v) 



A ,A . 

r’ r4-l 

Then rewriting (2.2) again 



(2.8) R = 21 fill ^ (u,T)du(lvdx 

rsrO 0 ^U 6 -x 4 rV ^1 *- V - U ^ 



Now, the problem is to find the values that give the ninlnun risk. 



Using classical optimization tecl'inique with Lagrangian fora the problem is to 



n 

Minimise ^.-1) 

n 

Subject to 21 ' 

J-0 J 
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I 



Since 



using X for the Lagrange mxxltiplier . 



(2.9)' $ = Z UJ [ z - Z ^ • A_1 f (u,v)dudvdx 

TirO O-U-Z-V-1 jrO V - U -i 



-^(Z A. - 1) 
0^0 ^ 



The partial derivative of ^ with respect to A is 



a| 




k-1 X - u \ z - u 

X A - 






J*-0 T- — U J V — U 



fj^(u,v)dudvdz 



- Z1 \ \ - Z A -A If (u,v)dudvdx +— = 0 

rrk>4-1 ^ ^ L j~0 ^ ■ V - uJ 2 

for k=1, , , , ,n-1i 

The partial derivative of f with respect to A and using (2.4) is 






J J y Y 



dvdx 



“fff rrl 3C-U 

Z_ \ \ \{3E - Z A -A )f (u,v)dudvdx 

r=l J -J i»0 ^ ^ V - u ^ 



+ Z=o 



The partial derivative of | with respect to ^ and using (2.5) is 
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aA 




n-l ^ x-u“\x-u . 

-2LA -A u”-' 

'* 1 - U J 1 - U 



A 

dudx f — - 0 
2n 



Let 



e(u,v) 



- u f r 

— . 1(1 ) = 12 f.(«» 



V - u 



r * ' r 

x^v-l 



v)dudv 



) i 2 f i Z 

Ij.(g) = g(u,v)f^(u,v)dudv , I^(g ) = g (u»v)f^Cu,v)dudv 



Using these notations , then the system of derivatives is 



ai 1 

(2.10) = / xL(g)dx 

0 ^ 



k-1 1 



1 



J Iv-Cg )dx ^ 

j=^ ^ 0 ^ ^ 0 ^ r=k+1 o’ 






xl (l)dx 



n r-1 1 n 1 \ 

H ZI A / I(l)dx- Z1 A_ / I(g)dx4-— =0 
r-ktl 0-0 0 r=k+1 0 2 



for k - 1 , , 



, n-1 



All terms in the above eguations reduce to Beta fxuictions . The relationship 
between Beta and Gamma ftmction 



1 . . r(m)r (n) (m-1)l(n-l)l 

B(m,n) = / v““\l - v)""^dv = = 

0 P (m+n) (nw-n-l)I 



is used heavily, V/e proceed to evaluate these terms . 
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Evaluation of 



J' xlj^(g)dx = 




nl 



dudvdz 



5 5 
V -u 



(k-.l)I(n-k-1)I 



k-1 , . »n-k-1 , , 

u (1 - v) dudv 



V — u 



4 ^ 2 2 

t. f r V -u 



k/. .n-k>-l, , 

u (1 - v) dudv 



2 1 1 V - u 



nl 



1 



Ck>1)l(n-k-1)l \ 3 



t 1 



k 0 



k+2/. xn-k-1, 

V Cl - v) dv 



1 1 



_ f _ Tf-'^'dv + — i - T)"-'‘-'dTl 



kf1 0 



1 1 

- f 

bl'2 0 



1 

2 



1 1,0 

— f v^'^^Cl 

kf1 0 



v)“-'=^’dT + — ' 



I T'^^Ct-vf-'^'dT 



k^2 0 



( 2 . 11 ) 



31w4 



6(n+l)(n+2) 
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Bvalrmtion of 



! 

f xl (1)dx- 1 1 1 dudvdx 

0 ^ J J J Cr-1 ) 1 (n-r-1 ) ! 



-- 1 - ^ \ I / 2 2v r-1,. sn-r-1, , 

i \ \ (v -u )u (1-v) dudv 



Cr-1 ) I Cn-r-1 ) ! 2 



(r-1 ) ! (n-r-1 ) ! 2 (r 0 






r4-2 0 



( 2 . 12 ) 



r+1 



(n+l)(n+2) 



Evaluation of 






nl 



(k-1 ) ! (n-k-1 ) ! 



dudvdx 




(k-1 ) ! (n-k-1 ) ! ( 2 -J J v-u 



k-1/, > n-k-1, , 

u (1-v) dudv 




/(1-v) 



n-k-1 



dudv 
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n! 



k+1 / . \n-k-1 , 

V (1-v) dv 



(k-1 ) ! (n-k-1 ) I 




kvi 0 J 0 



1 

(2.15) 

2(nt1 ) 



Evaluation of 

1 




(g )ix = 



J 




u 1 . , . n ! 

— y u (l-v)^~ ** dudvdx 

u/ (k-1 )! (n-k-1 ) I 




n! 

(k-1 ) ! (n-k-1 ) ! 







1 1 



J v‘=+'(t - V) 



k+1 0 



n-k-1 



dv 



] 



(2.14) = — ^ 

3(n4-1) 
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Evaluation of 






n! 



(n-r) ! (n-r-1 ) ! 



dudvdx 



n! 



1 .1 






(n-r) 1 (n-r-1 ) ! I r 0 



1 1 

( r-hl/. »n-r-1, 
J V (1-v) dv 

r+1 0 



(2.15) 



1 

n-»-1 



Evaluation of 



1 1 

J I^(e)dx is sane as evaluation of J Ij^(g)dr 

0 0 



(2.16) 



1 



j Ij,(g)<i:c = 



1 



2(n^-1) 



The use of these evaluations (2,11), (2,12), (2,15), (2.14), (2.15), 
(2.16) in (2,10), then the derivative with respect to is 



(2.17) .1^ 

6(n.2) 2 j-0 3 



1 n 

+ Z1 (r+1) 



(n-»-2) r-k-t-1 



n r-1 

Z ZA. 

r-k+1 j*0 



1 n X (nf 1 ) 

-- Z 

2 r=-k^-1 2 
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5k+4 



- 4 3(n4l) - 



5(k+2)(k+1) n 



- 6 



(n+2) r=k4l 



- 3 (1 - 4 3A(n4l) - 0 

or after evaluation the derivative (2.10) is 



c _ c 4 6 ^ C ^ 3X(n4l) 4 3n 
^ r^k+1 ^ 



3k +6k42 
n42 



for; k - 1 , . . . , n-1 



The derivative with respect to ^ is 
(2.19) 9^ 



9A. 



-JI 



X X 

(x ^) — n^l-v)^^ dvdx 

V V 



*i. 

r-1 js.0 ^ 



X - u V 

)f (u,v)dudvdx 4 — - 0 

0-y '' ^ V - u ^ 2 



•=-ttj d_v)“"^ dvdx - j (^)^ (1-v)“"^dvdx 



i (^yi) dx - E C^[l^(1)dx - Za fyg)dx 4~= 0 
rs-1 ' r=-1 f r»1 ' 



This can be reduced to 

n 2 

(2.2Q) - C. 4 6 Z C, = 3>(n4l) 4 3n -~ 

’ k=l n 42 



for ; k - 0 



The derivative with respect to y is 

(2.21) 9$ 



9 $ f r r X - u 

— = ir^. / A. .A — 
JJL -i.u 



X - u 



1 - u 



n-*1 ^ 

u dudx 4 - 0 
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can be reduced to 



af 2 

= A 5A(n+1) •= 0 

dA H4-2 ^ 



(2.22) C - 5^(n^-l) + for; k=n . 

nv2 

so in general 

^ 3k^+6k+2 

(2.23) Sc " ^ ^ ^r " + 5n — 

r=rk+l nv2 



for ; k s 0, . . 

Using (2.25) , then system of eguations in matrix form is 



5666 ... 60 








e - 2/(n 4 2) 


1566. ..So 




=2 




e - 11/(n 4 2) 


0 1 5 6 ... 6 0 

• • e • • • 




=3 

• 




e - 26/(n 4 2) 

• 


• • • • • 
• • c • 

0 0 ... 1 5 6 0 




• 

• 

Vi 




• 

• 

e - ( 3 n^- 6 n 42 )/(n 42 ) 


a 0 ... 0 1 5 0 




c 

n 




e - (3n^-1 )/(n+2) 


0 0 ... 0 0 1 -1 




_ Vl ^ 


'k 


6 - (3n^46n42)/(n42) 



Trtiere 6 = 5A(n+1) 4 5n 



n 
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After successive row subractions the matris is 



4 1 0 0 ... 0 0 








9/(n+2) 


1 4 1 0 ... 0 0 




=2 




15/(n+2) 


0 1 4 1 ... 0 0 

• • # # * 




s 

• 




21/(n+2) 

• 


• • • • 




• 


- 


• 


• • • 

0 0.. .1410 




• 

j=n.1 




• 

(6n-3)/(n42) 


0. 0 ... 0 1 4 1 




c 

n 




,(6n+3)/(n+2) 


0 0 . . . 0 0 1 -1 


; 


.Vi. 


1 


;€-(3n^4-6n-(.2)/(n+2)^ 



To solve this matrix e<^uation is considered the linear second-oi’der 



difference efiniation such that 



(2.24) 






6k 9 



n-+ 2 



with initial conditions C. - 0 and C . - 1 

0 n+1 



The corresponding homogeneous equation of (2.24) is 



(2.25) ‘ ♦ “k.2 = « 
heis the auxiliary eqiiation 

2 

n 4- 4m 4- t ss 0 

The roots of this equation can be calculated using the quadratic formul 
These roots are A - -2 + iT? and B - -2 - \T5*^ 
so the general sotion of (2.25) is given by 

(2.26) Cj^ = 4 + C* 
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I 



Trial solution of C, is 

k 



Using (2.24) 









1 



0 2(n+2) 






^ 1 



1 n+2 



then 



Q* _ 1 4 2k 

^ 2(n+2) 



So the general solution is 
(2.27) 



Ck - K2B^ +. 2(nf2) 



using initial conditions to determine the values of and ^2 



Cq » 0 in (2.27) implies 



(2.28) 






1 



2(im-2) 



=1 in (2.27) implies 



(2.29) 



K^A 



n4-1 






1 



2(m-2) 



Simxiltaneous solution of (2.28) and (2.29) sxe 






1 + A 



n^l 



K. •= 



' ■ 2(n+2) (a"' - b"*” 



and K_ - - 



2(n.2)(A”*-' - b”^' ) 



SO 



(2.30) 



= -_L_ ( k + ^ , (A^ .. B^) I 

•' (nr2) [ ^ 3 



for; k = 0, . . . , n+1 
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then 



using 



( 2 . 31 ) 




C 



k 




1 

(n+2) 



(1 - A)(a““^ 4 A^) - (1 - B)(b“"^ + 

f ox f k^O^«..«|n 
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.m.CALCULATICai OF RISK. 



( 5 . 1 ) 



Using, 



Prom eq,uation (2.8) the risk equation is 



Z1 W x^I (l)dx - 2C / xl (1) dx + cj /l (1)dz 
r=^0[0 0 0 



^6 ^ 



- 2A^ xl^(g)dx + 2C^A.^ f I^(g)dx 



I 

i 



x^I^(l)dx = 


(r+ 2 )(r+-l) 




r\ r» 


~(n4-5)(n+2Kn4-l ) 


X — 


V f • • • 9 Zl 


xl^d) dx * 


r-v -1 


r - 


- 0 1 • • • 9 ^ 


(ni-l)(n+ 2 ) 






1^(1 )dx 


1 

n -*-1 


r = 


Of • • • 9 ^ 


xl^(g)dx = 


5r + 4 


X ^ 


0 Tl 


6 (n+ 1 )(n+ 2 ) 




\J f • • efit 


Ij.(g)dx = 


1 






2 (n+l) 


r = 




w 9 • • • 1 n 


I^(g^)dx = 


1 


r = 


Of • • *911 


3(n^l) 
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So 



(in-1 j 



R 



3im4 



n 



rsO 



n 



1 

r=0 



"r+1 



C 

r 



n 



Calculation of Za; 



r-0 



Using (2.31) 



Z^r = — i-2 Z 5 1 „.| ' [d-A)(A°-^ + a'') . (1-B)(b“-''+ E^IT 

S) (nt2)^r.o [ 2 (a“"^' - b“*')L -JJ 



and using 
n 



^ (A^“^ + A^)= 2j 



A =2 



r=0 



rrO 



n 



1 - A 



1 - A 



n+1 



± t B--) = al s' 

rrO r-0 






= 2 



1 - B 



similiarly 
n 



(A 



2(n-r) _ ^2r^ = 2 21 A^ = 2 



2(n^.l) 



r=0 



1 - A 
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2 (n 4 -l) 



^ \ ^ 1 - B 

,2(n-r) . t,2t 'C- t,2t 



z ) = 2 Z = 2 ^ 

r -0 3 ^sO 1 - B 



and 



1 _ 1 - A 



1 - A"^ 1 - 6 (H-B) 



1 _ 1 - B 



1 - B rr? 6(1 4 A) 



1 - A _ 1 - B 



1 4 A 1 4 B 



0 _ = JL. (B°^^ - a“+^ 

1 + A 1 4 B 



then 

( 5 . 2 ) 



Za^ 



2 n-t -3 . ■ B 4-1 



?;0 (n+2)^ " (n+2)^ 



(1 -A)^A^ ^ (1 - B)^ - 12 

2 (a“'^^ - b"^+^ ) ^ 



-fT 



(a“+^ 4 b“*-^ 4 2 



_/ ^v2/ ,n^l ^n*^1 \ 

2 (n 4 . 2 ) (A - B ) 



n 



Calculation of Z 



r =0 



Using ( 2 . 30 ) 



1 n 



Z r*' 



c . c - 



rwO r^O 



f \ (r^ 4 2 r +-|) 4 (r f-i)| (n 42 )C^^^ 

r-5 I ^ 
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•V 



(r ^-| ) [(n+2) - r -ij 



1 ♦ b^3 [a”-^ - 4. 

4 (^n+1 _ jn+1j2 

here 

X (r^ + 2r 4--) - n(p4-1 ) (4n+U) -f 9(n4-1) 

r=0 ^ 12 



and 



“ • [a“'-''.a“ 



LA 

r=0 



- 4 B=^] [ A°-^ - a”' - 4 b"'*' ] 



(a“*-’ - b”*' >' 



24(a”**^*') 



■|^_2n43 ^ j2ii 43 _ (6n44)(A“'*'^ 4 b”'*'^) - (A^*' 4 B^*' ) 



-(6n47)(A” 4 b“> * ® <“*’>] 



Then 



n 



(3.5) .C_ » - 



n n 



r;rO 



r+1 r 



n(n+1 )(4n+14) 
12(n+2)‘ 



iJ(iill4J_L4X4 2lrC^rIcl 

n+2 L r«1 . ^ fr1 ^-1 



[a2”-\b^"* 5 -(6n44)(A”*2^-^).(A2».VB^.')_(6„+7)(A%unj^a(nH 

24(n42)^(A“*' 2 b“*'')^ 

Using (5.2) and (3.5) then 

2n^43n-1 f? (a”*' 4 E ) 

- 22(„42f (A-'-b5^ 

' ^ ;42)^(a-lb-v <“’> 

^.2n(A^^^ 4 C^'*'^) 4 8(:i4l ) 4 )*j 
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I?»C0MPAiU3GNS AIJD SIMIARY 



Prom (5.17) page 14 of Prof . ^.ead's paper 

L A 1 If 

(n+1) ^R(C) = J (n-1 )x(1-r)dz -f 2(n4.1 ) j \ (l-2)h(z) 
0 0 *- 



where 



A 

C = 



and 



n+1 



h(x) 



X / X - u \ 

= / — ) 

0 ' 1 - u ' 



2_ ^ ^ 

(n+l) k(C) + 4(n4l) I (1-x)h(x)dz 

6 0 



and evaluation of 



1 1 X/X-U V 

j (1-x)h(r)dx = / (1-z) /( ^jududx 

0 0 0 Vl-u./ 

f I ci-.)(-iu, 

0^u«i-x^l 1-u 



let 



then 



I4TX 



y = 



l-u 



0^y<^ 



X s u + (1 - u)y , 1-x=(l-u)(1— y) and dy = 



so 

A 1 

(4.1) (n+1)H(C) ~ 4 

6(n+.1) 3(n+2)(n4l) 

where 

R(C) is Pyke risk 



4 xh(l-x) 



udu 



dx 



1 - u 
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B jgV of sample c^^m^llative distribution function R(SDF) is 



1 I 1 

E(SDF) Si— I xCl-2)<ix * — 
no 6n 



then 

• n+1 

C4.2) (n+l)R(SDF) = 

,6n 

Optimum risk R(OPT) is derived in this paper is 



t4.3) Cn4l)RC0PT) = 






12Cn+2)‘ 



12Cn42)V’^-B^'-^ 



(2nv3) (A“4^P)+8(n-*-l ) )+2n(A^‘*'^+B“^^)+40(n4l ) 



24(n+2)^CA“^^-]^^^ 



In 4hSa eguation last two terms may be omitted. Since they are close to 
zero for n — 6 
so 

2n^<.3n-l 

(nvl )R(0PT) 5 — 

12(n+2)^ 



Cxuaputer output is ehoTna on page 25 for n=1, . . .,60. Equations 

(4.1) , (4.2) and (4.3) are used , 
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N 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 

60 



(N+1 )R(C) 



0.05555555 
0.08333331 
0.09999996 
0 . 11111100 
0. 11904750 
0. 12499990 
0. 12962950 
0. 13333320 
0. 13636360 
0. 13888870 
0. 14102550 
0.14285700 
0. 14444440 
0. 14583320 
0. 14705870 
0. 14814800 
0. 14912270 
0. 14999980 
0. 15079360 
0.15151500 
0. 15217380 
0.15277770 
0.15333320 
0.15384600 
0.15432090 
0.15476170 
0. 15517230 
0. 15555540 
0. 15591390 
0. 15624980 
0. 15656550 
0.15686260 
0. 15714280 
0. 15740730 
0.15765750 
0. 15789460 
0. 15811950 
0. 15833320 
0. 15853650 
0. 15873000 
0.15891460 
0. 15909080 
0.15925920 
0.15942010 
0. 15957430 
0.15972210 
0. 15986380 
0.15999990 
0. 16013060 
0. 16025630 
0. 16037730 
0.16049370 
0. 16060590 
0. 16071420 
0. 16081860 
0. 16091940 
0. 16101680 
0. 16111090 
0. 16120200 
0. 16129020 



(N+1)R(SCF) 



0.33333330 
0.25000000 
0.22222220 
0.20833330 
0.19999990 
0. 19444440 
0. 19047610 
0.18750000 
0.18518510 
0.18333320 
0.18181810 
0.18055 550 
0.17948710 
0.17857140 
0.17777770 
0. 177C8330 
0.17647050 
0.17592590 
0. 17543850 
0.17499990 
0. 17460310 
0. 17424230 
0.17391300 
0.17361110 
0. 17333320 
0.17307680 
0.17283940 
0. 17261900 
0. 17241370 
0.17222210 
0.17204290 
0.17187500 
0.17171710 
0.17156850 
0.17142850 
0.17129620 
0.17117110 
0.17105250 
0. 17094010 
0.17083320 
0. 17073160 
0.17063480 
0.17054250 
0. 17045450 
0. 17037030 
0. 17028980 
0. 17021270 
0.17013880 
0. 17006800 
0.16999990 
0.16993460 
0.16987170 
0. 16981130 
0.16975300 
0.16969690 
0. 16964280 
0.16959050 
0.16954020 
0. 16949140 
0.16944440 



(N+1 ) R(OPT) 



0.05353007 
0.07628661 
0.09249049 
0. 10353520 
0.11179110 
0. 11814010 
0. 12318110 
0.12727660 
0. 13066930 
0. 13352540 
0. 13596240 
0. 13806620 
0.13990060 
0. 14151430 
0. 14294460 
0. 14422110 
0. 14536740 
0.14640240 
0. 14734160 
0. 14819760 
0.14898100 
0. 14970070 
0.15036420 
0.15097770 
0. 15154670 
0.15207600 
0. 15256950 
0.15303060 
0. 15346260 
0.15386800 
0. 15424930 
0.15460860 
0. 15494760 
0.15526810 
0.15557160 
0. 15585930 
0.15613250 
0. 15639220 
0. 15663940 
0.15687500 
0.15709970 
0. 1 5731440 
0. 15751970 
0. 15771620 
0.15790440 
0.15808480 
0. 15825800 
0. 15842430 
0.15858420 
0.15873790 
0.15888600 
0. 15902860 
0.15916600 
0. 15929870 
0.15942670 
0.15955030 
0. 15966980 
0. 15978530 
0. 15989710 
0. 16000530 



LAST TWO TERMS 
OF 

(N+1 )R(OPT) 



-0.01649304 
-0.00857831 
-0.00582393 
-0.00399830 
-0.00294773 
-0.00225484 
-0.00178204 
-0.00144335 
-0.00119288 
-0.00100234 
-0.00085407 
-0.00073642 
-0.00064150 
-0. 00056382 
-0. 00049944 
-0.00044549 
-0.00039983 
-0.00036084 
-0.00032730 
-0.00029822 
-0.00027285 
-0.00025059 
-0.00023094 
-0.00021352 
-0.00019799 
-0.00018410 
-0.00017163 
-0.00016038 
-0.00015020 
-0.00014095 
-0.00013254 
-0.00012486 
-0.00011783 
-0.00011137 
-0.00010543 
-0.00009996 
-0.00009490 
-0.00009021 
-0.00008586 
-0.00008182 
-0.00007806 
-0.00007455 
-0.00007128 
- 0.00006821 
-0.00006534 
-0.00006265 
-0.00006012 
-0.00005773 
-0.00005549 
-0.00005338 
-0.00005138 
-0.00004950 
-0.00004771 
-0.00004603 
-0. 00004443 
-0.00004291 
-0.00004146 
-0.00004009 
-0.00003879 
-0.00003755 
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After conputer calculations for n-l,...»60 C4.1) t (4.2) and 
C4.5) is shown in (.4.4) and (4.5) 



(n^-1 )R 
.55 ■ 




* 

50 

(4.4) 

26 



10 



20 




C4.5) 



n 
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